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Chapter 10. Operators of the scalar Klein Gordon field 



10.1 The scalar charge-current density 

The Klein Gordon equation is generally associated with a scalar field, a 
field which values do not change under Lorentz transform, but only changes 
coordinates. We must realize that this is a special case however. Often 
more complex fields will have components which individually all obey the 
Klein Gordon equation, for instance the Dirac field which is a bi-spinor. 
The relevance of the Klein Gordon equation goes well beyond the scalar 
field with which it is generally associated. We can express a lorentz trans- 
form of the scalar field as. 



Where A is a 4x4 matrix representing the Lorentz transform. Note that 
A~^, the reversed transform, is used here. If the original field is concen- 
trated let's say at a position Xorg, then this coordinate will be transformed 
by A to a new position Xnew Working in the new reference frame and 
looking at Xnew we should obtain the value given by the original function 
tpixorg)- Therefor we need to reverse transform our coordinate Xnew back 



If we start with the particle's field in its rest frame and we transform to a 
frame in which it moves with a velocity v, then the volume containing the 
field will Lorentz contract by a factor 7. This means that the probability 
densities in this volume should increase by a factor 7 also. We have however 
defined our field defined as a scalar field which does not changes its values 
under transform. This means that can not represent the probability 
density of the scalar field as it does in case of the Schrodinger equation. 

We need a function which at one hand is proportional to 'tp*'ip, but at the 
other hand also increases with gamma. For the Schrodinger equation we 
have the expressions for the probability flux which actually do transform 
correct. If we denote the probability with P" and the probability flux 
with P* then we have a relativistically correct transforming 4-vector. The 
probability flux for the Schrodinger field is given by. 



(10.1) 



to X, 



org- 




(10.2) 
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This expression produces a real values per definition and it is instructive 
to substitute ip with a complex plane wave eigen-function. 

ip = ex.p(^-iEt/h+ip- x/h^ (10.3) 



We see that only the phase-change rate of tp ends up in the probability 
flux, while the imaginary part, representing the magnitude change of ip is 
eliminated by the subtraction. 

We can now see that we get the required expression for the probability 
density P° if we take the derivatives in time x°=ct. 

V" = ^(rp^-T^) (10.4) 



The operator (the timc-dcrivativc in this case) is applied symmetrically to 
both ip and ip* and the sign is adopted to select the imaginary part of the 
exponent representing the phase change rate in time. 

Since the total probability integrated over space is 1, and the total charge 
should be equal to e, we can use the above expressions as the charge- 
current density of the scalar Klein Gordon field in the Pauli-Weisskopf 
interpretation. 

J'^ - '^(r^P-'^A (10.5, 



2m \ dx^i dx/j, / 

The Pauli-Weisskopf interpretation easies the 'negative probability' issues 
associated with anti-particles given by. 

= exp (^ + iEt/h+ ip- x/h^ (10.6) 



In this case the sign of the charge-density changes indicating the opposite 
charge of the anti-particle. 

For a single particle wave-function we expect unitarity. This leads to the 
requirement that the four divergence of the charge-current density van- 
ishes. We need the continuity relation. 
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I ^ ^ _(n^*)^ + ^*n^_^^ =0 (10.7) 



^x^^ 2m \ dxt^ dx^ dx^ dxi^ 



That is. The Total current leaving an infinitesimal volume-element minus 

the current entering it gives the decrease of charge per unit of time in the 
volume element. If obeys the Klein Gordon equation then we can replace 
this expression by: 



— ' mipip + ipmip— — — - — I = (10.8) 



^x^^ 2m \ dxi^ dx^ dx/j, dx^^ 



Which trivially vanishes. This means that J'* is a conserved current and 
that unitarity is guaranteed. 



10.2 Klein Gordon charge-current density operators 

In the previous section wc did sec that we can define the electric charge- 
density operator and the current-density operator for the Klein Gordon 
field as follows. 

~, ^ iehj_^ (10.9) 
2m dx° ' 2m dx^ 



Recall that j° = cp as being part of a four-vector. For this we need to 
define the way these operators act on the field as. 

j° = = ( ip* y ip + y* ip* ) 

^ ; ( (10.10) 

J' = -i/;* y t/; = - ( ^* J' V + ^ J'* V'* ) 

The operators act symmetrically on both ip and The practical reason 
is that the individual terms can contain imaginary parts and we expect 
the end result to be real. The imaginary parts cancel if we let the operator 
act in the above symmetric way. 



10.3 General form for the local operator application 
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The imaginary terms occur as soon as ip is not a plane wave solution 
anymore, when its magnitude changes with time or space, for instance 
because it is localized. 

The imaginary parts also vanish if we integrate the terms over space. 
3m!^J dx^ii;* r i;^^ = 3m dx^ (^^^ r*tp*^ ^ = 0(10.11) 



They cancel because any variation in the magnitude of V can be seen 
as the result of interference between plane waves which do have constant 
magnitudes them self. Any interference between plane waves is a sinusoidal 
function which vanishes after being integrated over space. 

The fact that we want to evaluate quantities locally instead of globally is 
what leads us to the expressions of (10.10) 



10.3 General form for the local operator application 

The charge-current density operators suggest that a more general form for 
operators should be. 

O = ^**0 V = ( 4^*0 V' + ^ O*^* ) (10.12) 

If the field ip is locally described by an exponent with an arbitrary ampli- 
tude part a{x'^) and an arbitrary phase <p{x'^) like this: 

(p = exp(^a(x^) +i4>{x'')^ (10.13) 

Then the general operator form of (10.12) guarantees that the change in 
amplitude doesn't play a role. Only the phase change rates of ip do. There 
is a good reason to assume that for an important class of operators only the 
phase change rates determine the quantities associated with the operators. 

For instance if we consider a static, typically bound, solution at rest with- 
out internal momentum then the phase change rate in time determines 
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the energy. Changing reference frames will transform this energy into mo- 
mentum via the Lorcntz transform. The class of local quantities for which 
we want to have the operators are typically based on the relativistic four- 
momentum. 

The local quantities arc densities and generally this is where the local 
amplitude of the wave-function comes in, in the form of the conjugate 
product adapted by a factor 7 due to the Lorentz contraction. The 
latter is specific for the scalar Klein Gordon field because the value of cf) 
itself, being a Lorentz scalar, does not transform. 

We will find that (10.12) is not the most general form. The problem arises 
when the operator contains higher order/ mixed derivatives. In these cases 
the phase and amplitude components do mix up again and we have to study 
the results carefully to find the correct way to apply the operators in these 
cases. This is the price we pay for representing a field with both phase and 
amplitude by complex numbers. 



10.4 The time derivative Hamiltonian 

The use of the time derivative Hamiltonian is often used to derive new 
operators from existing operators. More specifically, the operators which 
give the time derivatives of the quantities produced by existing operators. 



The commutation of the operator with the Hamiltonian is used to obtain 
the time derivative of the observable. 



Although taking the time derivative directly is trivial, applying H is gen- 
erally much more interesting since it is an alternative representation which 
should lead to the same results. Often H doesn't contain any time- 
derivative at all but obtains the same result by other means. 




(10.14) 




(10.15) 



10.4 The time derivative Hamiltonian 
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Naively applying [O, H] however leads easily to erroneous expressions, es- 
pecially when we want them to be locally meaningful rather then globally 
after taking the integral over all space. We will therefor review here the 
conditions under which we are allowed to use the commutator with the 
Hamiltonian. 

The time derivative of the local density of a quantity is given by. 

|(r6*) = ^-olW + + §^{r)di, (10.16) 



The more general operators do not change with time them shelf. This 
eliminates the middle term at the right hand side. Substituting the time 
derivatives with H gives us. 



1 

ih 



(10.17) 



Where we have used H* = —H in the second occurrence of H at the right 
hand side since it is Hermitic. We see that the correctness of the use of 
the commutator of the Hamiltonian as a time derivative depends on the 
presumption of. 

? 

[H*iP*^dil) = iP*Hd'^ (10.18) 



This assumption however can be made only in a limited number of cases, 
even if the operator O is unity! The magnitude of the field I'^l has to be 
static for the expression below to be valid. ( = ) 

{H*ip*)ij = iIj*HiIj (10.19) 



The above identity is not true if the time derivative of tjj isn't purely 
imaginary. This condition is relaxed however if we integrate over all space. 

J dx^ H*tp*ij = Jdx^ij*Hij (10.20) 
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This identity is valid in general for linear H. In this case 4> may be an 
arbitrary superposition of plane waves. The sinusoidal interference terms 
between the plane waves integrate away. 

= Jdx^ ^tVi + V'2V'2 + V'sV's + -) (10.21) 



The complex conjugate product of the superposition is equal to the sum 
of the conjugate products of the individual components. 

tp = exp(^a(a;'') +i^/)(x^)) (10.22) 

10.5 The position operator 

For a position operator we can simply define a center of gravity operator. 
Starting from the charge-density operator, removing the charge from the 
normalization and multiplying with the position 

2mc ox° 

The operator is applied in the general form. 

= V* ? V' = ( V'* X* V + V' X" V* ) (10.24) 



(10.23) 



The density is compressed by a factor 7 (higher by a factor 7) due to 
Lorentz contraction. It is the time derivative applied anti-symmetrically 
on -0 and ■0* which produces this factor 7. The phase change rate in time 
is proportional to the energy E = ^mc?. 

= a;^7V'*V' (10.25) 



The integral over space gives the average position, the center of gravity. 

Jdx^tp*¥ij = J dx^ x'-fip*!^ (10.26) 



•^avg 



10.6 The velocity operator 
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10.6 The velocity operator 

We would like, as in non-relativistic quantum mechanics, to obtain the 
velocity operator from the position operator via a commutation with the 
Hamiltonian. 

(10.27) 




However the Klein Gordon equation doesn't explicitly contain the first 
derivative in time but the second order instead which gives us the square 
of the Hamiltonian. 

= -^^^V' = ( -nVV^V + ?"^c^ JV' (10.28) 



One way to obtain the Hamiltonian is via a derivation against itself. 

1 d 

H = -^(^ -h^c^V'^^p + n?c^^^ (10.29) 

and using the time derivative ihdt for H at the righthand side. The ratio- 
nale for this is that many operators, like for instance the position operator 
depend on the density which involves a time derivative. 

The only non-commuting part of the Hamiltonian with the position oper- 
ator is the second order derivative of the corresponding axis, so. 



1 

ih 



ih d .2 2^^ 



2mc dx° ' dx"^ 



h^c d d 
m dx dx" 



The time derivative in X was handled just like a constant here. We ob- 
tain the velocity operator after applying the derivative with regard to the 
Hamiltonian. 

V- = l^f5!f (10.31) 

2 dH \ m dx dx° J ^ ' 

As said, to do this we define the Hamiltonian H as the time derivative 
here. This then removes the time derivative from the expression and we 
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obtain the velocity operator. Which turns out to be proportional to the 
current-density operator in (10.9). 

~i ih d , , 

V = -7^1^ 10.32 



Applying the operator on the Klein Gordon field in the standard way, 
V' = V* V* V = ( V'* V + ^ v"v* ) (10.33) 

gives us F*, the velocity-density of the Klein Gordon field 

^ 2m i ^ dx^ dx^ ^ / (10.34) 



If we compare this with (10.5), for the current density, we see that they 
are equal up to the value e for the electric charge. 

/ = (10.35) 



Which justifies the term velocity density for V*. The integral of over all 
space gives us the velocity of the particle, or more accurate, the spatial 
average velocity. 

When wc apply the velocity operator on the Klein Gordon field tp expressed 
in a form which separates amplitude and phase. 

ijj = exp{a + i(t)) (10.36) 



We obtain the following expression for the velocity density. 

V' = -14^*^ (10-37) 
m ox^ 



Where the only dependence on the amplitude component a is found in 
The subtraction in (10.34) has eliminated the da/dx terms, so the velocity 



10.6 The velocity operator 
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density docs not depend on the derivatives of a, which is indeed what we 
want. The spatial phase change rate is proportional to the momentum. 

d(f> p jmv' n^Qe^ 

So we can express the velocity density V* in a form which uses the local 
velocity and the corresponding local 7 

= jv' (10.39) 



We have defined the Klein Gordon field V' as a Lorentz scalar and thus 
the term does no change under Lorentz transform. The total volume 
which contains the field however does change. It undergoes a Lorentz 
contraction by a factor 7. This means that we have the following expression 
for the unitarity of the field. 

J dx^ = 1 (10.40) 



Where the local 7 compensates the local Lorentz contraction. We can now 
show that the integral over all space indeed leads to the (spatial average) 
velocity of the field. 

V'^^g = J r dx^ = J v'jtlj*ipdx^ (10.41) 



or expressed using (10.34) 



We have established a way to obtain a local velocity f*, with a corre- 
sponding local 7, from the Klein Gordon field using the velocity density 
operator V*. 

The average velocity is obtained by averaging the local velocity 
over space using the relativistic effective density 7^*^. 
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The local 7 can be obtained from the three components of 7?;* via the 
standard relativistic equation. 



7 = cosh(sinh \p\) = cosh ( sinh ^ J {'jVx)'^ + (jVyY + {iv^Y ) (10.44) 



Nevertheless we need to be aware that what we call the local v and local 7 
are in fact average values in case of superposition and the more the wave- 
packet is localized the higher is the amount of superposition and the wider 
the range of momenta. 

The local velocity we obtain is the same velocity we obtain if we use the 
current density divided by the charge density. Classically, in the case of 
a current carrying electric wire, this is called the average drift velocity. 
In an actual wire there are many charges moving at all kinds of different 
speeds. 

We can however carry out the electromagnetic calculations for the Lorentz 
force and transform them correctly from one to another reference frames 
only with the knowledge of these averaged current and the charge, without 

the specific knowledge of the speed of all the individual charges. (Sec the 
appendix on the Magnetism as a relativistic side effect of electrostatics 
where the mechanism behind this is explained) 



10.7 The acceleration operator 

Having derived the velocity operator, leading to a standard result, we like 
to go one step further here: We want to obtain the acceleration operator. 
We do so by applying the Hamiltonian twice on the position operator. 



10.7 The acceleration operator 
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The rightmost form is tlic result of tlic cancelation of two oppositely signed 
terms of the form HXn. To work out the above we use the square of the 
Hamiltonian operator. 



9t2 



( -n^c^VV + m^c^^il) (10.46) 



The mass term m?c^ commutes with the position operator, so it does 
not contribute. The non-commuting terms we need are the second order 
derivatives over the axis corresponding with the ?-th component of the 
position operator X*. 



1 



ih 



d 



X 



2„2 



2mc dx° 



-hfc 



d_ 



ihc d d 
m dx° dx^ 



(10.47) 



When we compare this with (10.32) for the velocity operator V* we see 
an extra time derivative. This is what we would expect from the classical 
case. 



Meaning of the acceleration operator 

If we apply the operator A* in the standard way on the Klein Gordon field, 
then we might hope to get an acceleration-density A* in the form of. 

? 

= a'j^*4) (10.48) 



Where a* is the local acceleration, and the factor '^ip*%l} is the local rela- 
tivistic density. Since is a Lorentz scalar we need a local 7 here to 
correct the density for the local Lorentz contraction. Compare this to the 
position-density and velocity-density we obtained. 

= x' r = v'jiIj*i/j (10.49) 



Since f '7 is proportional to the the momentum, the phase change rate over 
space, and we take the time-derivative: 
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d_ 
di 



V 7 



( 7^ a' jip*^; ) (10.50) 



We do not get the average acceleration if we integrate A* over all space. 
So what do we get? Recall the definition of the spatial components of the 
relativistic 4-velocity 



It is the velocity measured in the proper time of the particle. The factor 
7 is valid for an arbitrary accelerating particle since. 

T = f^dt =^ ?r = ^ (10-52) 



7 dt J 



The relativistic four-acceleration is the second derivative of x* in regard 
with the proper time r. 

^ = -wTjr: = -WT- 10.53 

OT Ot Ot Ot J 



Wc now can see that A* is the product of the relativistic density 7 
the relativistic four-acceleration. 

A' = ^7V'*V' (10.54) 



Which means that the integration over all space of A' gives us the average 
relativistic four-acceleration. 



5V 



avg 



J A* dx^ = J dx^ (10.55) 
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Application of the acceleration operator 

If we express the Klein Gordon field in a way which separates the phase 
and its amplitude, and furthermore define related fields ^ and |V'P as. 

■0 = exp {a + i4>) 

= exp(z(/)) (10.56) 
IV-P = exp (2a) 



Then we can apply the operator O^, operating on the phase only, in the 
following way. 

O' = (^O'* (10.57) 

As long as contains no higher as first order derivatives then we can 
write this in the usual asymmetric form. 

IV'P = - (10.58) 



This won't work however if the operator contains higher order derivatives 
as is the case of the acceleration operator. When we apply the operator A* 
in the usual anti-symmetric way on the Klein Gordon field ip, 

7^ --(r^^>P - ^^^r) (10.59) 



then we get the amplitude component a mixed in the result. The applica- 
tion of A* that correctly does what we want is given by evaluating (10.50). 



d 
d 



ih 

2m 



d f 1 dip 1 dip* 



dt \ ip dx^ ip* dx^ 



ip*ip (10.60) 



ih f ^^d d ^ _ ^ d d _ V'* d^j ^ i) d^f dip* 
2m \ dt 9x* dt dx^ ip dt dx'- ip* dt dx'- 
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This more cumbersome expression has two additional terms which cancel 
the unwanted terms which are depending on the amplitude. It can be 
checked, by substituting ^ = exp(a + icp) that this expression is equiva- 
lent to, 

as in equation (10.57), where = exp(i^). Our confidence in (10.60) is 
boosted considerably because it can be used to derive the acceleration from 
the electro magnetic fields, (as will be shown briefiy in the next section). 
Working out (10.61) gives us. 

dt dx^ dt ^ ^ ^ (^^) 

The result does not depend on the order in which the derivatives are ap- 
plied, but note that on exchanging dx^ and dt we get. 

m ( dx^ dt^ ^ ^ m ( dx^ ^ V' "0 (10.63) 

Because we associate d(p/dx^ with the momentum but we associate 
d(f)/dt with the energy p° instead. The phase (f)\s s, scalar function of the 
four component which leads to a redundancies. 

The above shows that there are two terms which contribute identically to 
(f) and thus to A' also. We have to consider them both. Since the two 
expressions are equal we can write. 

¥ = (10.64) 
dt dx' ^ ' 



This tells us that there will be an acceleration in the direction of x' if the 
(potential) energy decreases along x*, which can for instance be the case 
in classical general relativity near a massive body. 

In this case the field is an interacting field and the potential difference 
should actually be described by an interaction term with gravitational 
field. 



10.8 Acceleration operator with EM interaction 
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10.8 Acceleration operator with EM interaction 

The chapter on the " Lorentz force derived from Klein Gordon's equation" 
handles the electromagnetic interaction with the Klein Gordon field in 
great detail. Here we will briefly preview how the 4-0'(^(^^l^''^o,t'>'On density 
operator we have just arrived leads to the electromagnetic Lorentz force. 

The four-potential A'^ determines the electromagnetic interaction. Mul- 
tiple like charges brought close to each other increase each others energy 
(phase change rate), while multiple charges moving in a bunch increase 
each other's momentum (spatial phase change rate) . 

The total phase change rates determines the total four-momentum, the so 
called canonical momentum. This is the sum of the inertial momentum 
(from the mass) plus the interaction momentum. The physical velocity is 
still determined by the inertial momentum. The effect of the electromag- 
netic four-potential is an additional phase change added to the field when 
propagating. 

If A^^ is zero and we increase it then it doesn't simply adds phase change 
rates to the field. Rather the field will change it's inertial momentum to 
compensate the added phase change rates. The simple reason is special 
relativity. The phase of the field at a certain point would depend on 
the spatial integral of the phase rates. If A^^ changes at an arbitrary far 
away point then the phase would change anywhere instantaneously. Hence, 
inertial momentum has to change in order to compensate the induced phase 
change rates from the four-potential. 

The derivatives determine the total canonical momentum. We use the 
gauge covariant derivative D^, which removes the electromagnetic interac- 
tion to recover the inertial momentum which determines the actual physical 
speed 



We therefor replace the normal derivatives in our expression for the four- 
acceleration density (10.60) with the covariant ones. 



D,j, = dn + ieA^/h 



(10.65) 




d 



d 



e 




d 



d 



e 



(10.66) 
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Working this out step by step gives us the following expression. 



ox" ax* 
dx° ax* 
^ dx° dx' ^ 



{i — 



ihc 
2m 



- i-ib A^—— 
dx° 



^ ^ 9x* 

' fi^ 9x* 

+ 'ft"* 



(10.67) 



nr 



+ -^A°A'ij*^ 



The first column here is simply (10.60) without the interaction. The fourth 
column adds up to zero and the third column adds up to zero as well. 
The second column contains the difference between interaction and no 
interaction. 



If we express the total canonical momentum by capitol then the inter- 
action version of (10.62) becomes simply. 

dP^ ( dp' dA'\ 



It doesn't make any difference if we change the order of the two derivatives. 
The interacting version of expression (10.63) becomes. 

e'^) (10.70) 



dx^ \ dx^ dx^ 



The two expressions are degenerate with respect to each other since is a 
scalar function determined by two vectors and which determine the 
phase (f). The represent two causes with the same outcome and we have to 
consider them both. 
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Equating the two gives us the expression for the four- acceleration due to 
electromagnetic interaction. 

M = _ ^ , (_^_^'\ nn7n 

dx'' dxi V dx" dxi ^ ' ' 



The last term contains the definition of the Electric field E* so that we can 
write. 

= - 1^ + ^ P (10.72) 

ox° ox^ c 



We have found the Electric part of the Lorentz force. We obtain the 
Magnetic part if we consider the phase changes at a moving point instead 
of a point at rest. This is handled in great detail and with lots of illustrating 
images in the chapter " Lorentz force derived from Klein Gordon's equation" 



10.9 The angular momentum operator 

The angular momentum operator should extract the quantum mechani- 
cal equivalent of the classical f x p. To get the density of the angular 
momentum we need the density of the momentum. 

More to come here... 



